Optical one-way quantum computing with a simulated valence-bond solid 
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One-way quantum computation proceeds by sequentially measuring individual spins (qubits) in 
an entangled many-spin resource state [1]. It remains a challenge, however, to efficiently produce 
such resource states. Is it possible to reduce the task of generating these states to simply cooling a 
quantum many-body system to its ground state? Cluster states, the canonical resource for one-way 
quantum computing, do not naturally occur as ground states of physical systems [ , ]. This led 
to a significant effort to identify alternative resource states that appear as ground states in spin 
lattices [4-8]. An appealing candidate is a valence-bond-solid state described by Affieck, Kennedy, 
Lieb, and Tasaki (AKLT) [ ]. It is the unique, gapped ground state for a two-body Hamiltonian 
on a spin-1 chain, and can be used as a resource for one-way quantum computing [ - ]. Here, we 
experimentally generate a photonic AKLT state and use it to implement single-qubit quantum logic 
gates. 



In the standard circuit model of quantum computa- 
tion [in], information is carried by two-level systems 
called qubits. The computation proceeds dynamically 
via unitary single-qubit logic gates and multiple-qubit 
entangling gates. Apart from these entangling gates the 
qubits are fully isolated from each other. Computations 
in the one-way model, on the other hand, are performed 
via single-qubit measurements on a strongly-correlated, 
i.e., entangled, resource state. The one-way model has 
led to some of the highest estimated error thresholds for 
fault-tolerant quantum computation [ , jlz], and to a 
series of experimental demonstrations of quantum logic 
gates [13-18], wherein the technical requirements can be 
much simpler than for the circuit model. This is particu- 
larly true of optical implementations, where the resource 
requirements for one-way quantum computing are sig- 
nificantly lower [ ], and the predicted error thresholds 
significantly higher [ ], than for any other approach to 
quantum computation. 

Because qubits in the one-way model are not isolated 
but rather interact strongly with each other, this ap- 
proach lends itself more naturally for implementations 
in condensed-matter systems. But, out of the vast va- 
riety of strongly-coupled quantum many-body systems, 
can we find one that has a ground state we can use as a 
resource for quantum computing? That seems unlikely if 
this ground state is to be the cluster state, because the 
cluster state is not the ground state of a strongly-coupled 
many-body system with a Hamiltonian consisting of two- 
body interactions [ , ]. As a result, the search for al- 
ternative resource states has attracted a lot of interest 
recently. Although up to now little is known about the 
requirements potential resource states for the one-way 
model have to meet, and although most states are in fact 
useless for this task [ ], a handful of alternative states 
have been identified [4-5]. All of these states can be 
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FIG. 1: AKLT states, (a) The AKLT state [ ] is a valence- 
bond solid and can be represented by a chain of spin- 1 singlet 
states where adjoining qubits of neighbouring pairs are pro- 
jected on the triplet subspace, i.e. the subspace symmetric 
with respect to swapping of the two qubits. At either end of 
the chain a boundary qubit remains, ensuring that the ground 
state is non-degenerate, (b) One can simulate an AKLT state 
with a chain of sources producing singlet states and projecting 
pairs of particles on the triplet subspace. 



described in the framework of projected entangled pair 
states [ , ] or matrix product states [ , ] . 

A promising candidate is the ground state of a spin 
model studied by Affleck, Kennedy, Lieb, and Tasaki 
(AKLT) [ ]. This valence-bond-solid state (see Fig- 
ure la) appears as the unique gapped ground state 
of a rotationally-invariant nearest-neighbour two-body 
Hamiltonian on a spin-1 chain. The AKLT state pos- 
sesses diverging localisable entanglement length [ >] and, 
remarkably, can serve as a resource for one-way quantum 
computation [6, 7, 24]. Because the Hamiltonian is frus- 
tration free, i.e. the ground state minimises the energy of 
each local term of the Hamiltonian, measurements in the 
course of the computation leave the remaining particles in 
their ground state. Operations leaving the computational 
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FIG. 2: Tomographic reconstruction of our photonic 
AKLT state, (a) and (b) show the real and imaginary part 
of the density matrix reconstructed from an over-complete 
set of tomography measurements. The fidelity with the ideal 
AKLT state is (87.1 =b 0.4)%. 



subspace are penalised by the energy gap protecting the 
AKLT state. Universal quantum computation can be 
achieved via dynamical coupling of several AKLT states, 
where each can be regarded as 'quantum computational 
wires' [o, Y, 24]. These properties render the AKLT state 
an attractive alternative to cluster states as a more natu- 
ral resource for quantum computing in condensed-matter 
systems. 

Quantum computation with AKLT states is differ- 
ent from computing with cluster states in a number of 
ways. The elementary physical units are spin-1 systems 
(qutrits) instead of spin-^ systems (qubits), although it 
is still qubits that are encoded as 'logical' information. 
Adaptive measurements allow the performance of non- 
Pauli operations, including Clifford gates. Single-qubit 
rotations can be performed around any Cartesian axis. 
These operations are probabilistic, rather than determin- 
istic, and succeed with probability |. When an operation 



fails, it performs a heralded logical-identity operation, 
i.e. a teleportation of the logical information along the 
chain. The operation can then be reattempted on the 
next qutrit until it succeeds. Combinations of such rota- 
tions allow the implementation of arbitrary single-qubit 
quantum logic gates. 

Although a number of one-dimensional spin chains are 
well-described by the AKLT Hamiltonian, most promi- 
nently Ni(C2H8N2)2N02(C104) (NENP) [ ], up to now 
experimental techniques do not allow the single-spin mea- 
surements necessary for one-way quantum computation. 
Yet, one of the fundamental and most appealing mo- 
tivations for quantum computing, is the possibility to 
simulate aspects of quantum systems that cannot di- 
rectly be studied [26]. Because the AKLT state is a 
valence-bond solid state (see Figure la), we can sim- 
ulate it via a chain of spin-^ singlet states, for exam- 
ple polarisation-entangled photon pairs, where adjoining 
particles of neighbouring pairs are projected on the sym- 
metric triplet subspace (see Figure lb). While this ap- 
proach does not allow to analyse the dynamics of the 
corresponding solid-state system, it allows the direct pro- 
duction of the AKLT state, and to use it for one-way 
quantum computation. 

Here, we experimentally demonstrate the generation of 
photonic AKLT states and their application for one-way 
quantum computation. We produce two singlet states, 
= (|HV) — |VH))/>/2, in four distinct spatial modes 
using spontaneous parametric down conversion. Here, 
|H) and |V) denote horizontal and vertical polarisation, 
our computational basis. From these two singlets we cre- 
ate an AKLT state consisting of two boundary spin- 1/2 
systems and one spin-1 system. The spin-1 system is a 
biphoton, symmetrised by projecting a pair of photons 
into the triplet subspace. Qutrit measurements are per- 
formed using the method from Ref. [ ~]. For details of 
our experimental setup, see the Methods section and Ap- 
pendix A. A detailed discussion of the theoretical aspects 
of the simulation of AKLT states using quantum optics 
and their use in one-way quantum computation will be 
given in a separate paper [28]. 

The AKLT state for a qubit-qutrit-qubit system is 
\i^AKLT) = ^|H,1,V) + ^|V,1,H)-^|H,2,H)- 
^|V,0,V). In our case, the qutrit states |0), 
1 2) correspond to the biphoton states ^<^H<^H|vac), 
ttfja^lvac), -^a\^a\^\vdic) ^ respectively. Here, |vac) is the 

vacuum state, and and are photon creation oper- 
ators. To verify the faithful production of the AKLT 
state in our experiment, we perform quantum-state to- 
mography on the qubit-qutrit-qubit system and use a 
maximum-likelihood technique based on a semi-definite- 
programming algorithm [29] to reconstruct the density 
matrix shown in Figure 2. For a detailed list of the states 
measured and of the corresponding counts, see the Table 
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FIG. 3: Measurement results for single-qubit rota- 
tions.. (a)-(c) show the coordinates of the Bloch vectors of 
the reconstructed output density matrices for rotations of a 
logical input state |H) around the x, y and z axes, respec- 
tively. Note that the results shown are for the 'plus' outcome 
of the qutrit measurement, and that we have applied the nec- 
essary Pauli corrections to the reconstructed density matrices 
for all plots shown in the figure. Error bars are standard 
deviations calculated from Monte-Carlo simulations. Solid 
and dashed lines indicate the theoretical expectations given 
the ideal AKLT state and the tomographically-reconstruced 
AKLT state (see Figure 2), respectively. For the rotation an- 
gles 0, ^, and |, panels (d)-(f) show the Bloch vectors of the 
measured (and Pauli corrected) density matrices correspond- 
ing to the Bloch coordinates shown in (a)-(c). 



V in Appendix B. The fidelity of the reconstructed den- 
sity matrix with the ideal AKLT state is (87.1 ± 0.4)%. 
We calculate this value in a Monte-Carlo simulation with 
420 iterations on the observed counts, using the defini- 

/ / \ 2 

tion F(p, cr) = Tr a/ ^/ap^/a ] for the fidelity between 

two quantum states [ ]. 

In order to demonstrate the use of AKLT states for 
quantum computation, we realise single-qubit rotations 
around the f , y and z axis of the Bloch sphere. We begin 



the computation by projecting the first boundary qubit 
onto some qubit state By doing so, we effectively 
prepare the logical state IV^^), where (V^IV^^) = 0. To 
perform a rotation Rx{0)^ ^yi^) o^" Rz{^) of the logi- 
cal state by an angle around the respective coordinate 
axis, we project the qutrit on one of the corresponding 
bases given in Table L Qutrit measurements are realised 
probabilistically as described in the Method section and 
Appendix A. We denote the three outcomes of each qutrit 
basis as 'plus', 'minus' and 'id'. Each is expected to occur 
with probability 1/3. Up to a known Pauli error [ , , ], 
which can be corrected as indicated in Table I, the out- 
comes 'plus' and 'minus' signal a successful rotation, and 
the outcome 'id' signals the logical identity, i.e., a rota- 
tion by 0°. As a result, a successful rotation is achieved 
with probability 2/3. For ^ = 0, every outcome heralds 
the logical identity. This can be used to teleport logi- 
cal information along the wire, for example to a position 
where the wire is coupled to another, or to the read-out 
position. 

To prepare our logical input state, we project the first 
qubit on one of a set of states: |H), |V), |()^), 
|m±). Here, |±) = (|H) ± |V))/V^, are the eigen- 
states of the Hadamard operator, |m+) = cos(|)|H) + 
e*?sin(|)|V) is the 'magic state' [ ] with ^ = 
arccos(^), and (m~|m+) = 0. For each axis of rotation 

we choose 10 angles = {O, f , J, f , f , f f , f , ^} 
and project the qutrit on the corresponding state (see 
Table I) for the 'plus' and 'minus' outcomes. We project 
on the 'id' outcome once for every input state and ro- 
tation axis. Finally, we reconstruct the density matrix 
of the computational outcome by performing a tomo- 
graphically over-complete set of measurements on the 
last qubit, using the measurement settings: |H), |V), 
|±), |R) = (|H)+i|V))/V2, and |L) = (|H) - i|V))/V2. 

Figure 3 shows measurement results for single-qubit 
rotations of the logical input state |H) (i.e. projecting 
the first qubit on |V)) around the three rotation axes. 
The plots in Figure 3a-c) show the coordinates of the 
rotated Bloch vectors as compared with the theoretical 
expectation. In Table II we list the fidelities for rota- 
tions of |H) as well as the averaged fidelities for all logi- 
cal input states prepared. The probabilities for the three 
qutrit measurements, averaged over all input states, rota- 
tions and rotation angles, are measured to be 0.34 ±0.03, 
0.30±0.05, and 0.36±0.04 for the 'plus', the 'minus', and 
the 'id' outcome, respectively. This is in good agreement 
with the expected value of | for each outcome. An aver- 
age of the output fidelities achieved over all input states 
and all rotations performed yields a value of (92 ±4)%, 
demonstrating the high quality of our single-qubit quan- 
tum logic gates using a photonic AKLT state. A detailed 
list of all results can be found in Appendix B. 

We have experimentally demonstrated a one-way 
quantum-computation scheme harnessing a novel re- 
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TABLE I: Qutrit measurement bases and Pauli corrections. Single-qubit rotations are realised by a projective measure- 
ment in a corresponding qutrit basis that has three possible outcomes: 'plus', 'minus' and 'id'. The qutrit states \y)^ and 

1 2;) are defined as | (^Q^hCIh ~ '^v'^v) I vac), | (^Q^hCIh + ^v^v) |vac), and ^ ^^h^^v + '^v'^h) I vac), respectively, and X, y, Z 
indicate the Pauli correction that has to be applied to the read-out qubit depending on measurement outcome and measurement 
basis. 



gate fidelities for logical input H) 


outcomes 


Rx 


Ry 




Pth Pexp 


Pth Pexp 


Pth Pexp 


plus 


0.91 ±0.04 0.98 ±0.02 


0.90 ±0.05 0.98 ±0.02 


0.90 ±0.03 0.98 ±0.02 


minus 


0.93 ±0.03 0.97 ±0.03 


0.91 ±0.03 0.99 ±0.01 


0.92 ±0.04 0.97 ±0.02 


id 


0.90 ±0.03 0.98 ±0.02 


0.92 ±0.02 0.999 ±0.006 


0.97 ±0.02 0.99 ±0.02 




gate fidelities 


averaged over all input states 


outcomes 


Rx 


Ry 


Rz 


Pth Pexp 


Pth Pexp 


Pth Pexp 


all 


0.92 ±0.04 0.97 ±0.02 


0.91 ±0.04 0.98 ±0.01 


0.92 ±0.04 0.98 ±0.02 



TABLE IL Single-qubit logic gate fidelities. We compare the experimentally determined output density matrices with the 
ones expected, on the one hand, given an ideal AKLT state, pth^ and, on the other hand, given the AKLT state measured in our 
setup, pexp- The upper part of the table shows the fidelities for a logical input state |H). For the 'plus' and 'minus' outcomes 
the fidelities are averaged over all rotation angles, for the 'id' outcome we performed one measurement per rotation axis. The 
lower part shows the corresponding fidelities averaged over all logical input states prepared (see text) and over all three qutrit 
measurement outcomes. 



source, the AKLT state, and used it to implement a cir- 
cuit realising single-qubit rotations around any coordi- 
nate axis. Quantum computation using AKLT instead 
of cluster states promises to combine the inherent ad- 
vantages of the one-way model with resources that oc- 
cur naturally in physical systems. Our scheme for cre- 
ating AKLT states uses entangled states and linear op- 
tics similar in requirements to optical implementations 
using cluster states [19]. In contrast to some other opti- 
cal implementations of quantum logic gates for one-way 
quantum computation [ , ], our scheme does not re- 
quire phase stability and achieves significantly higher ex- 
perimental fidelities. Our implementation of a valence- 
bond-solid state is a realisation of a projected entangled 
pair state [22]. Such states offer a promising framework 
for understanding the properties of entangled states that 
make them useful computational resources [4-6, 8]. Gen- 
eralisations of the presented approach might allow to sim- 
ulate other classes of alternative resource states with lin- 
ear optics and to study their potential for quantum com- 
puting. Future challenges will be to find efficient meth- 



ods of coupling quantum wires, to study solid-state com- 
pounds with ground states that can be used as computa- 
tional resources, and to implement techniques to address 
such systems on a single-particle level. Ideally, this and 
related research will lead to implementations in solid- 
state architectures, allowing to tap the power of one- 
way quantum computation while taking full advantage 
of the appealing characteristics of novel resource states 
like AKLT. 

We thank W. A. Coish, R. Prevedel, A. C. Doherty and 
A. Gilchrist for valuable discussions. We are grateful for 
financial support from NSERC, OCE, CFI, Quantum- 
Works, and MRI ERA. 

Methods 

The light source in our experiment is a Tita- 
nium: Sapphire femtosecond laser, centred at 790 nm with 
10 nm full-width-at-half-maximum (FWHM) bandwidth, 
2.9 W average output power and 80 MHz repetition rate. 
Second-harmonic generation in a 2 mm thick Bismuth- 
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Borate crystal (BiBO) yields a beam of 780 mW power, 
centred at 395 nm, with about 1.5 nm FWHM band- 
width. With this beam we pump two separate type-I 
spontaneous parametric down- conversion sources [31, 32], 
each a pair of 1 mm thick beta-Barium-Borate (BBO) 
crystals. Longitudinal and transverse walk-off occurring 
in the down- conversion crystals is compensated with a 
combination of birefringent crystals (a-BBO, quartz and 
BiBO, see Ref. [ ] and Appendix A). All photons pass 
through 3nm FWHM bandwidth filters. In each source, 
the polarisation of the photons in one mode is measured 
directly at the source, the photons in the other modes are 
coupled into single-mode fibres and sent to a quantum in- 
terferometer and analyser setup. The input modes of the 
interferometer are overlapped at a 50:50 beam splitter 
(BS), where, depending on the two-photon state, two- 
photon interference leads to both photons leaving via the 
same or via different BS output modes [ ]. By mea- 
suring a two-photon event in one output mode of the 
BS, the biphoton is projected onto a qutrit subspace. 
This mode is input in a qutrit analyser, where we im- 
plement qutrit projections by probabilistically separating 
the two photons at another BS and performing appropri- 
ate polarisation measurements on each photon [27]. If we 
assume that the two polarisation measurements project 
on the states |V^rn) = ^cosa^a^j + e*^"^ sina^a^^ |vac) 

(m = 1,2), we can calculate the qutrit state this mea- 
surement projects on by propagating these states back 
through the BS. For a more detailed discussion of the 
setup and the qutrit projections, see Appendix A. 

Appendix A: Experimental Setup 

The entangled photon pairs in our experiment are gen- 
erated using two separate type-I spontaneous parametric 
down conversion (SPDC) sources [31, 32]. Each consists 
of a pair of 1 mm thick beta-Barium-Borate (/3-BaB204, 
or BBO) crystals, their optical axes oriented perpendic- 
ular to each other. Longitudinal walk off in the SPDC 
crystals is compensated using a 0.5mm quartz, a 2mm 
quartz and a 1mm a-BBO crystal before the first SPDC 
source, and a 2mm a-BBO and a 1mm quartz crystal 
between the two sources. Additional transverse walk-off 
is compensated by placing 1mm thick BiBsOe crystals 
cut at ^ = 152.6° and (/) = 0° in modes 2 and 3 (see Fig- 
ure 4a). These angles are chosen such that the crystals 
compensate transverse walk-off without introducing lon- 
gitudinal walk off. All photons pass through 3nm FWHM 
bandwidth filters. The phases in the setup and the polar- 
isation rotation in the single-mode fibres is set such that 
the sources produce singlet states in modes 1 & 2 and 3 & 
4. In modes 1 & 2 we measure a fidelity of (96.9 ± 0.5)% 



with the ideal singlet state = ^(|HV) - |VH)), 

and a tangle of 0.92 ± 0.01. For the second source the 
fidelity is (96.9 ±0.6)% and the tangle is 0.90 ±0.01. We 
had single count rates of about 200kIIz in the qubit anal- 
ysers 1 and 2 (see Figure 4a), and single count rates of 
around 80kHz in the detectors Di and D2 in the qutrit 
analyser (see Figure 4c) (both sources contribute to these 
latter single count rates). The two- fold coincidence count 
rate for the first source was 7.4kHz between qubit anal- 
yser 1 and Di in the qutrit analyser. For the second 
source we had a two- fold coincidence count rate of 5.9kHz 
between qubit analyser 2 and Di in the qutrit analyser. 

One photon of each pair is measured directly at the 
source, using polarisation analysers. The modes for both 
measurement outcomes are coupled into single-mode fi- 
bres and monitored via single-photon detectors (Perkin- 
Elmer, SPCM-4Q4C). The two remaining photons are 
coupled into single-mode fibres and sent to a quantum 
interferometer and analyser setup. The input modes of 
the interferometer are overlapped at a 50:50 beam split- 
ter (BS). If the input photons are set to have the same 
polarisation, Hong-Ou-Mandel (HOM) interference [34] 
occurs. Postselecting on four-fold events with one pho- 
ton in mode 1, one in mode 4 and two photons in the 
output mode of the BS indicated in Figure 4a, we ob- 
serve constructive HOM interference with a visibility of 
95.7 ±3.7%. 

Whether two photons entering the interferometer (via 
modes 2 and 3) leave through the same or different BS 
outputs, depends on the biphoton state [ ]. In partic- 
ular, coincidence detection events between two different 
BS outputs only occur for the two-photon singlet state. 
By post-selecting on a biphoton excitation in one output 
mode of the BS, the biphoton is projected onto a symmet- 
ric subspace and can be described as a qutrit [ ]. This 
is identical to the symmetrisation needed to generate an 
AKLT state [ ] (see Fig. 1 in the main text). We measure 
this qutrit using the analyser outlined in Figure 4c. 

The analyser works by probabilistically splitting up the 
biphoton at a BS, and by performing a qubit projective 
measurement on each of the output modes 1 and 2 of the 
BS. We project on a given qutrit state by projecting the 
two photons on corresponding qubit states. In particular, 
the photon in mode m = 1,2 of the qutrit analyser is 
projected on state IV^m) = cos a^|H)±e*^'^ sina^|V). A 
successful projective measurement of a qutrit is heralded 
by a coincidence event between detectors Di and D2. 
To calculate which qutrit state such a coincidence event 
signals we can propagate our two-qubit state |?/^2) 
back through the beam splitter. All contributions in the 
unused of the two input ports of the BS can be neglected. 
The (unnormalised) two-photon state then becomes: 
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FIG. 4: Experimental Setup, (a) Two type-I SPDC sources are used to generate entangled pairs. Longitudinal and transverse 
walk-off are compensated via combinations of birefringent crystals Tl, T2 and C. All photons are coupled into single-mode 
fibers. Polarisation rotation in the fibres is compensated via polarisation controllers, and the phase is adjusted by tilting quarter 
waveplates (QWP). Photons in modes 2 and 3 are fed into a 50:50 beam splitter (BS). Detecting two photons in one of the 
BS outputs projects these photons on the symmetric qutrit subspace. They can then be treated as a qutrit. (b) Each qubit 
analyser consists of a half waveplate (HWP), a QWP and a polarising BS. Both outputs of the polarising BS are monitored, 
(c) In the qutrit analyser the biphoton forming the qutrit is probabilistically split up using a BS. Phases introduced at the BS 
are compensated by tilted QWPs. A combination of QWPs and HWPs allows to project the qutrit onto a state of our choice. 
A successful projection is heralded by a coincidence event between detectors Di and D2. 
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TABLE III: Analyser parameters for qutrit measurements for rotation gates. In the qutrit analyser, one photon is projected on 
I -01), the second one on \ip2), where \iprn) = cosarn|H) + 6*^"^ sinarn|V). Each qutrit measurement has three possible outcomes, 
'plus', 'minus' and 'id' corresponding to three different sets of parameters for the qutrit analyser. We provide the settings for 
rotations around each of the Cartesian axes. 
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cos ai cos 0^2 a|ja|j + (e*^^ cos ^2 sin ai + e*^^ cos ai sin ^2) a^^a^ + e**^^^+^^^ sin ai sin a2a^a^ |vac). (AI) 

In our qutrit basis we can write this as: 

—zz cosai cosa2|0) H — (e*^^ cosa2 sinai + e*^^ cosai sin 0^2) |1) H — Q'^ixi+X2) ^[yiq:^ sina2|2). (A2) 
V 2 2 y 2 



r 



In general, the success probability of this projective 
qutrit measurement depends on the qutrit state. For ex- 
ample, to project on the biphoton state a|ja|j|vac), we 
set both analysers to |H). Given that the biphoton is in 
the correct state, a coincidence will occur with probabil- 
ity ^ because of the probabilistic splitting of the photons 
at the BS. As a second example, in order to project on 
the biphoton state a|ja^|vac), we choose {i/ji) = |H) and 
|'02) = |V). If the biphoton is in the right state, the suc- 
cess probability will only be | because the photons can 



be split in four possible ways, and only one leads to a 
coincidence event. 

Table III lists the parameters am and Xm we have to 
choose to project on a qutrit measurement for a rotation 
by an angle 6 around axis x, y or z. The probability for 
each of these qutrit measurements to work is |. By tak- 
ing into account the | V) outcomes of the polarising beam 
splitters (PBSs) in the qutrit analyser, this probability 
could be increased to |. In order to analyse the qubit- 
qutrit-qubit state generated in our setup we performed a 
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TABLE IV: Analyser parameters for qutrit measurements used for over-complete tomography. Each qutrit measurement is 
implemented by splitting the biphoton probabilistically at a beam splitter and projecting the two photons (m = 1,2) on 

iV^m) = cosam|H) + e*^"^ sin ttrn I V). For brevity, we use the definitions ^ = arccos(^) and 77 = arccos(Y^), and p denotes 
the success probability for a given qutrit projection. These success probabilities are taken into account by our tomographic 
reconstruction of the density matrix. 



tomographically over- complete set of measurements and 
performed tomography using a maximum likelihood tech- 
nique 35]. The set of measurements performed for 
the two qubits was |H), |V), |+), |-), |R), and |L). In 
Table IV we list the set of qutrit measurements as well 
as the corresponding parameters am and Xm for the set- 
tings of the qubit measurements in the qutrit analyser. 
The table also shows the probability with which each of 
the projections succeeds given that the biphoton is in the 
corresponding qutrit state. 

A qubit analyser where both outputs (|H) and |V)) 
are monitored can be operated in various ways. One is 
to adjust the coupling efficiencies of the two single- mode 
fibre couplers such that they are approximately the same. 
This method is, however, very susceptible to long-time 
drifts of the sources because these might again lead to 
an unbalance between the two couplers. The method we 
chose is the following. Let us assume we want to measure 
the counts for a projection on |H). Then we record the 
counts in both PBS outputs once with the waveplates set 
such that the transmitted path in the PBS corresponds 
to |H), and we record the counts in both PBS outputs 
for the same amount of time with the waveplates set such 
that the transmitted path corresponds to |V). We do the 
same for any polarisation we project on, and by adding up 
the respective counts we average over any imperfections 
in the balance of the coupling efficiencies. If we have N 
qubit analysers for which we apply this technique, the 
number of combinations to measure will be 2^. For our 
tomography measurements of the AKLT state we used 
this technique for modes 1 and 4. For the single-qubit 
rotation measurements we used it for mode 4, but we 
only monitored the transmitted mode in analyser 1 in 
order to simplify the software implementation of the scan 
routine. 



Appendix B: Results 

All measured counts correspond to four-fold coinci- 
dence detection events between one detector in each of 



the two qubit analysers and the two detectors Di and 
D2 in the qutrit analyser (see Fig. 4). Table V lists the 
four-fold coincidence counts measured for all tomograph- 
ical settings. For each setting of the analyser waveplates 
we integrated over 60s. Because we monitored both out- 
puts in each of the two qubit analysers and applied the 
technique described above in order to average over any 
unbalance between the two analyser outputs, we have to 
measure 4 combinations of settings per projective mea- 
surement. This results in 240s overall measurement time 
per projective measurement. In order to reduce the ef- 
fect of slow drifts in the setup, we performed the full set 
of measurements twice, in each case randomly ordering 
the settings, resulting in a measurement time of 480s per 
setting in Table V. 

In the main text we gave the measurement results for 
the rotation of a logical input state |H) around all three 
coordinate axes. For completeness, we have done the 
same measurements for a set of logical input states. It 
is important to note, that in our setup it is possible to 
rotate an arbitrary input state around any of the three 
coordinate axes because we project the first qubit on a 
given state rather than on a measurement basis. As we 
mentioned in the main text we prepare a logical input 
state by projecting on the qubit state orthogonal to it, 
i.e. IV^^). Because we post select on four- fold coincidence 
events, we automatically disregard those cases where the 
outcome of our measurement would be If one takes 
into account both outcomes of the projective measure- 
ment, each of them will occur randomly with probability 
^, and if we get the outcome |?/;), the logical input state 
will be the state orthogonal to what we want to prepare. 
To correct for this error for any arbitrary state is impos- 
sible because that would require a universal-NOT opera- 
tion, which is non-unitary [36]. In practice, one chooses 
input states such that errors in the preparation can be 
corrected via Pauli operations, i.e. any of the states along 
the coordinate axes of the Bloch sphere. 

For each input state, rotation angle, qutrit outcome 
and rotation angle we performed a tomographically over- 
complete set of measurements (|H), |V), |+), |— ), |R), 
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^/2 


V2 
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3 


4 
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141 
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62 


41 
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85 
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|V) 
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11 


13 


319 


55 


96 


80 


43 


156 
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70 


67 


70 


61 


1+) 


127 


173 


193 


18 


231 


5 


292 


117 


21 


138 


162 


144 


155 


69 


120 




|-> 
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34 


168 


34 


85 


84 


42 


6 


140 


79 


52 


52 


56 


71 


83 




|R) 


126 


85 


139 


23 
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51 


156 


74 


84 


52 


184 


21 


192 


14 


151 




|L> 


122 


88 


205 


27 


180 


41 


162 


75 


84 


160 


49 


166 


35 


149 
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|H> 


8 


64 


320 


80 


54 
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23 


67 


66 


52 


68 


149 
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66 
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315 


87 


15 


303 


16 


67 


43 


61 


99 
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147 
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57 


80 
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56 
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8 


286 
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14 
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75 
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87 


148 
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23 
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80 
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56 
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30 
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157 


11 




|L> 


154 


89 


167 


152 


26 


179 


32 


79 


92 


196 


39 


204 


45 


20 


125 




|H> 


10 


58 


270 


41 


42 
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139 


77 


84 


39 


45 


34 


300 


96 
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|V) 
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85 


5 


161 


119 


59 


71 
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68 


19 
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76 


65 


73 
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1+) 


116 


68 
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32 


140 


26 
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67 


99 


12 


153 


52 
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159 


93 
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57 
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34 


94 


71 


24 


192 


55 


168 
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8 




|R) 


135 


126 


183 


172 


128 


156 


149 


12 


137 


8 


255 


15 


296 


87 


89 




|L> 


116 


22 


192 


57 


46 


93 


68 


162 


18 


35 


66 


96 


62 


57 
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|H> 


9 


63 


229 


60 


37 


142 


134 


97 


66 


56 


53 
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24 


47 


71 




|V) 
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73 


6 


165 


183 


58 


67 


84 


86 
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20 


86 


45 


99 


63 


|L> 


1+) 


144 


85 


167 


36 
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36 


214 
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72 
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46 
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37 


143 


13 
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128 
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49 
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30 


97 


80 


221 


23 
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37 


20 
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|R) 
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13 
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61 


75 


57 


71 
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8 
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43 


45 


62 


68 


93 




|L> 


133 


152 


172 


149 


124 


150 


163 


10 


131 
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14 


306 


15 


82 


82 



TABLE V: Table of tomography results. We performed a tomographically-overcomplete set of measurements on our qubit- 
qutrit-qubit state. The states in the first column indicate the measurement settings for the first qubit, the state in the second 
column indicate the settings for the other qubit, and the states in the top row denote the qutrit measurement settings. All 
counts given are raw four-fold coincidence events integrated over 480s. We have not corrected these counts to take into account 
the probabilities of success for the various qutrit measurements involved. In order to make this correction to get the "real" 
count rates, which one can use for the tomographic reconstruction, one has to divide the raw counts by the probability of 
success for the qutrit measurement as given in Table IV. 



and |L)) on the qubit carrying the result of our single- 
qubit logic gate. In Tables VI, VII, and VIII we list the 
fidelities of the reconstructed density matrices for all ro- 
tations of various input states compared with what we 
would expect given a perfect AKLT state. Table VI 
shows these fidelities for rotations Rx{0) around the x 
axis by an angle 0. For each of the 6 different logical in- 
put states and each angle we observed the 'plus' as well as 
the 'minus' outcome of the qutrit measurement. Because 
the 'id' outcome is independent of the rotation angle, we 



only measured it once for every rotation axis. For each 
measurement we performed a Monte Carlo simulation on 
the measured counts, which are assumed to be the mean 
of a Poissonian count distribution, to reconstruct a set 
of 400 density matrices. These we used to calculate the 
means and standard deviations for the fidelities given in 
this table. Corresponding results for rotations Ry{0) are 
given in Table VII. For rotations Rz{0) around the z axis 
(see Table VIII) we performed measurements for two ad- 
ditional logical input states, 
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id n.a. 0.90 zb 0.02 0.90 zb 0.02 0.88 zb 0.04 0.96 zb 0.02 0.90 zb 0.03 0.93 ± 0.03 



TABLE VI: Fidelities for rotations around the x axis. This table lists the fidelities for various input states of the computation 
compared with the theoretical expectation given an ideal AKLT state. These fidelities are given for a number of angles of 
rotation around the x axis. We list them separately for the three distinct outcomes 'plus', 'minus' and 'id' of the qutrit 
measurement. For each input state we measured the 'id' outcome only once because it is independent of the rotation angle. 
The value and error for each fidelity are the mean and standard deviation from Monte Carlo simulations based on a Poissonian 
distribution around the counts measured. 
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zb 


0. 


,03 


0, 


.90 


zb 


0. 


,02 


0. 


,95 zb 


0. 


,02 


0. 


,97 


zb 


0, 


.02 


0. 


,95 


zb 


0, 


.02 


37r/2 


0, 


.87 


zb 


0, 


.03 


0, 


.88 


zb 


0, 


,02 


0, 


.87 


zb 


0, 


,04 


0, 


,93 zb 


0, 


,03 


0, 


,95 


zb 


0, 


.03 


0, 


,90 


zb 


0, 


.03 


77r/4 


0, 


.92 


zb 


0, 


.03 


0, 


.90 


zb 


0. 


,03 


0, 


.91 


zb 


0. 


,02 


0. 


,93 zb 


0. 


,02 


0. 


,93 


zb 


0, 


.03 


0. 


,99 


zb 


0, 


.01 





0, 


.90 


zb 


0, 


.02 


0, 


.86 


zb 


0. 


,02 


0, 


.89 


zb 


0. 


,04 


0. 


,86 zb 


0. 


,03 


0. 


,94 


zb 


0, 


.03 


0. 


,91 


zb 


0, 


.03 


7r/8 


0, 


.91 


zb 


0, 


.03 


0, 


.89 


zb 


0. 


,03 


0, 


.90 


zb 


0. 


,03 


0. 


,92 zb 


0. 


,02 


0. 


,98 


zb 


0, 


.02 


0. 


,98 


zb 


0, 


.02 


7r/4 


0, 


.93 


zb 





.03 





.90 


zb 


0, 


,03 





.92 


zb 


0, 


,02 


0, 


,93 zb 


0, 


,02 


0, 


,94 


zb 


0, 


,02 


0, 


,90 


zb 


0, 


,03 


37r/8 


0, 


.91 


± 


0, 


.03 


0, 


.89 


± 


0, 


,02 





.92 


± 


0, 


,03 


0, 


,89 ± 


0, 


,03 


0, 


,91 


± 


0, 


,03 


0, 


,95 


± 


0, 


,02 


7r/2 


0, 


.91 


zb 


0, 


.02 


0, 


.91 


zb 


0, 


,02 


0, 


.96 


zb 


0, 


,03 


0, 


,85 ± 


0, 


,03 


0, 


,91 


zb 


0, 


.03 


0, 


,96 


zb 


0, 


.02 


37r/4 


0, 


.90 


zb 


0, 


.04 


0, 


.96 


zb 


0. 


,03 


0, 


.93 


zb 


0. 


,02 


0. 


,89 zb 


0. 


,03 


0. 


,95 


zb 


0, 


.02 


0. 


,93 


zb 


0, 


.02 


TT 


0, 


.91 


zb 


0, 


.02 


0, 


.94 


zb 


0. 


,02 


0, 


.89 


zb 


0. 


,03 


0. 


,92 zb 


0. 


,04 


0. 


,95 


zb 


0, 


.02 


0. 


,95 


zb 


0, 


.02 


57r/4 


0, 


.95 


zb 


0, 


.02 


0, 


.99 


zb 


0, 


,01 


0, 


.89 


zb 


0, 


,02 


0, 


,84 zb 


0, 


,03 


0, 


,91 


zb 


0, 


.03 


0, 


,96 


zb 


0, 


.02 


37r/2 


0, 


.92 


zb 


0, 


.02 


0, 


.86 


zb 


0. 


,02 


0, 


.84 


zb 


0. 


,04 


0. 


,82 zb 


0. 


,03 


0. 


,90 


zb 


0, 


.03 


0. 


,91 


zb 


0, 


.03 


77r/4 


0, 


.83 


zb 


0, 


.04 


0, 


.90 


zb 


0. 


,03 


0, 


.92 


zb 


0. 


,02 


0. 


,91 zb 


0. 


,02 


0. 


,86 


zb 


0, 


.03 


0. 


,95 


zb 


0, 


.02 



id n.a. 0.92 zb 0.02 0.88 zb 0.02 0.89 zb 0.04 0.93 zb 0.03 0.95 zb 0.03 0.93 zb 0.03 



TABLE VII: Fidelities for rotations around the y axis. 
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logical input state 












|H> 


IV) 




ir> 


|m+) 


|m-> 


1+) 


|-> 



plus 



0.87 zb 0.03 0.94 zb 0.02 0.89 zb 0.04 0.89 zb 0.04 0.93 zb 0.03 0.91 zb 0.03 0.92 zb 0.03 0.90 zb 0.03 

7r/8 0.92 zb 0.02 0.91 zb 0.02 0.86 zb 0.04 0.92 zb 0.03 0.92 zb 0.03 0.93 zb 0.03 0.97 zb 0.02 0.98 zb 0.01 

7r/4 0.85 zb 0.03 0.85 zb 0.03 0.91 zb 0.04 0.89 zb 0.04 0.94 zb 0.03 0.94 zb 0.03 0.97 zb 0.03 0.89 zb 0.04 

37r/8 0.85 ± 0.03 0.95 ± 0.01 0.88 zb 0.04 0.89 dz 0.03 0.86 zb 0.04 0.91 zb 0.03 0.92 ± 0.03 0.88 ± 0.04 

7r/2 0.90 ± 0.03 0.92 ± 0.02 0.79 ± 0.06 0.98 ± 0.01 0.85 ± 0.04 0.86 ± 0.03 0.88 ± 0.05 0.89 ± 0.04 

37r/4 0.86 zb 0.04 0.93 zb 0.02 0.95 zb 0.04 0.93 zb 0.04 0.83 zb 0.04 0.88 zb 0.04 0.92 zb 0.04 0.87 zb 0.06 

TT 0.92 zb 0.03 0.87 zb 0.03 0.79 zb 0.05 0.84 zb 0.05 0.91 zb 0.04 0.90 zb 0.04 0.95 zb 0.03 0.96 zb 0.02 

57r/4 0.91 zb 0.03 0.87 zb 0.04 0.98 zb 0.01 0.91 zb 0.05 0.89 zb 0.03 0.99 zb 0.01 0.93 zb 0.04 0.96 zb 0.04 

37r/2 0.90 zb 0.04 0.95 zb 0.02 0.92 zb 0.05 0.90 zb 0.05 0.94 zb 0.03 0.88 zb 0.04 0.91 zb 0.03 0.93 ± 0.03 

77r/4 0.94 zb 0.03 0.99 zb 0.01 0.91 zb 0.06 0.99 zb 0.01 0.93 zb 0.03 0.96 zb 0.03 0.83 zb 0.06 0.96 zb 0.04 



0.87 zb 0.03 0.95 zb 0.02 0.88 zb 0.04 0.88 zb 0.03 0.90 zb 0.03 0.98 zb 0.02 0.90 zb 0.03 0.92 zb 0.03 

7r/8 0.89 zb 0.02 0.90 zb 0.03 0.93 zb 0.03 0.96 zb 0.03 0.94 zb 0.03 0.94 zb 0.03 0.92 zb 0.03 0.88 zb 0.04 

7r/4 0.95 zb 0.01 0.90 zb 0.02 0.92 ± 0.04 0.89 dz 0.04 0.96 zb 0.02 0.97 zb 0.03 0.89 zb 0.05 0.94 zb 0.04 

37r/8 0.91 ± 0.02 0.93 ± 0.02 0.96 ± 0.03 0.92 ± 0.04 0.90 ± 0.03 0.94 ± 0.03 0.90 ± 0.04 0.94 ± 0.04 

7r/2 0.90 zb 0.03 0.89 zb 0.03 0.96 ± 0.03 0.85 zb 0.06 0.90 ± 0.04 0.91 zb 0.03 0.91 zb 0.03 0.94 zb 0.03 

37r/4 0.93 zb 0.02 0.94 zb 0.02 0.99 zb 0.02 0.93 zb 0.04 0.91 zb 0.04 0.98 zb 0.01 0.97 zb 0.03 0.99 zb 0.01 

TT 0.86 zb 0.04 0.87 zb 0.03 0.93 zb 0.04 0.92 zb 0.04 0.90 zb 0.04 0.90 zb 0.04 0.92 zb 0.03 0.90 zb 0.04 

57r/4 0.88 zb 0.05 0.99 zb 0.01 0.94 zb 0.05 0.89 zb 0.05 0.89 zb 0.04 0.94 zb 0.03 0.90 zb 0.05 0.80 zb 0.06 

37r/2 0.96 zb 0.03 0.91 zb 0.04 0.83 zb 0.09 0.98 zb 0.02 0.79 zb 0.04 0.86 zb 0.03 0.97 zb 0.02 0.96 zb 0.02 

77r/4 0.98 zb 0.01 0.99 zb 0.01 0.93 zb 0.05 0.87 zb 0.06 0.92 zb 0.03 0.83 zb 0.04 0.88 zb 0.05 0.94 zb 0.04 



id 



n.a. 0.96 zb 0.01 0.97 zb 0.01 0.96 zb 0.03 0.96 zb 0.02 0.93 zb 0.03 0.98 zb 0.02 0.90 zb 0.03 0.90 zb 0.03 



TABLE VIII: Fidelities for rotations around the z axis. 



